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We have performed first order reversal curve measurements of the temperature-driven metal- 
insulator transition in VO2 thin films, which enable quantitative analysis of the hysteresis behavior. 
An unexpected tail-like feature in the contour plot of the reversal curve distribution indicates the 
existence of metallic domains, even at temperatures below the closing of the hysteresis. These 
domains interact with the surrounding medium and change the reversal path relative to a path from 
a fully insulating state. With this in mind, and assuming that such interaction persist through the 
entire phase transition, we develop a model where the driving force (or energy barrier) in charge of 
opening a hysteresis in VO2 are inter-domain interactions. This model is intrinsically different from 
the Preisach model usually used to describe hysteresis; given that it looks for the microscopic origin 
of the hysteresis, and provides physical parameters to characterize it. 

PACS numbers: 71.30.+h, 72.80.Ga, 64.60.Ht,64.60.an 



I. INTRODUCTION 

Vanadium oxide (VO2) exhibits a sharp first-order 
metal-insulator transition (MIT) and structural-phase 
transition (SPT) from a high temperature metallic rutile 
phase to a low temperature insulating monoclinic phaso^, 
with a hysteresis of ~1K in bulk samples^ii^ which can 
reach a few degrees in thin films^i^. 

The nature of the electronic and structural phase tran- 
sition has been a matter of research and debate for over 
three decades, with many recent experiments shedding 
additional light on the issue^ii^i^ii^iiiii^iHiiiiis^ But 

during this period the hysteretic nature of VO2 has re- 
ceived much less attention. 

The hysteresis nature and mechanism in VO2 is still 
an open question which is interesting from a basic re- 
search view point, and might provide both information 
regarding the MIT mechanism in VO2 and of hysteresis 
in general. It is also of vast importance from a tech- 
nological aspect. For example, for optical memory-type 
applications, a large hysteresis is needed, while a small 
hysteresis is preferable in VO2 based room temperature 
bolometers for UV detectors. So, understanding the ori- 
gin of the hysteresis can lead to a better control of its 
properties. The VO2 thermal hysteresis characteristics 
are found to vary between samples grown under differ- 
ent conditions, such as thickness of thin filmsiS,, growth 
temperatureii, grain sizei^, and choice of substrate^J^. 
For example, in polycrystalline thin films and single crys- 
tals of VO2 embedded in silicon, the hysteresis width 
grows with decreasing of the crystallite sizoi^. Attempts 
in modeling the thermal hysteresis have been mainly by 
various forms of the Preisach model^S, where hysteresis 
is assumed to exist in different domains of the system un- 
der investigation. Here the hysteresis widths and ampli- 
tudes, and distribution of domains are fitting parameters 
of the modelSSi'^. To the best of our knowledge, there is 



a scarcity of physical models of the hysteresis and a dis- 
cussion of its microscopic origin, especially for first order 
phase transitions. 

A powerful tool for studying the hysteresis in var- 
ious system is the first-order reversal curve (FORC) 
analysis^Sii^. FORC diagrams are based on the pro- 
cedure first-described by MayergoyaSi. They are used 
extensively in the study of magnetic systems, providing 
information regarding different phases^^, interactions^^ 
and magnetic state of the system, examples include the 
identification of a transition from single domain to a vor- 
tex state in nanoscale magnets2L2^, and one dimensional 
bubble states which exhibit large coercive fields in Co/Pt 
multilayers^^. The FORC method can be extended to 
any system where the response to an external parameter 
shows hysteretic behavior—. This has been utilized in the 
study of ferroelectric materials^, spin-crossover solids^ 
and spin glasses^^. 

In this paper we report on our FORC measurements 
of the temperature driven MIT in VO2 thin films. We 
prove the FORC analysis to be a useful tool in quan- 
tifying the hysteresis of different samples, which is im- 
portant in any attempt to correlate the hysteresis pa- 
rameters with various sample characteristics. Moreover, 
we find that even at low temperatures (down to 32 IK), 
where the VO2 films seem to have no hysteresis, there 
is an irreversible signal in the FORC analysis. This fea- 
ture is attributed to persistence of metallic domains down 
to these low temperatures, and indicates that there are 
interactions between the metallic domains and the sur- 
rounding films during the phase transition. Assuming 
that such interactions should exist throughout the phase 
separated transition of the VO2 filmiii^, we propose a 
phenomenological model of the hysteresis in VO2 sam- 
ples. We take the interactions between the different VO2 
domains as the driving force behind the opening of the 
hysteresis, which is fundamentally different from the var- 
ious forms of the Preisach model, since it does not assume 
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the existence of hysteresis in every domain. This model 
enables us to reconstruct the entire FORC diagram from 
one full heating branch and two fitting parameters. 

In Sec. HI] the experimental setup and procedures are 
described in detail. In Sec. Illll we show the experimental 
data, along with a physical interpretation. In Sec. II VI we 
describe our theoretical model of thermal hysteresis, and 
Sec. |V]is devoted to a summary of our findings. 



II. EXPERIMENTAL 
A. Sample preparation 

Vanadium oxide thin films were prepared by reactive 
RF magnetron sputtering of a vanadium target (1.5" di- 
ameter, 99.8%) on r-cut (1012) sapphire substrates. The 
samples were prepared in a high- vacuum deposition sys- 
tem with a base pressure of 5xl0~^ Torr. A mixture of 
ultra-high purity (UHP) Argon and UHP oxygen gasses 
were used for sputtering. The total pressure during de- 
position was 3x10"'^ Torr, and the oxygen partial pres- 
sure was optimized to 1.5x10^^ Torr (5% of the total 
pressure). The substrate temperature during deposition 
was 500 °C; while the effective RF- magnetron power was 
kept at 300W. These conditions yielded a deposition rate 
of 0.37 A/s. The samples were cooled at a rate of 13 
°C/min in the same Ar/02 flow of the deposition. Films 
were characterized and verified to be single phase VO2 by 
X-ray diffraction using CuKq, radiation and by Energy- 
dispersive X-ray spectroscopy. Surface morphology was 
measured by Atomic force Microscopy (AFM) and Scan- 
ning Electron Microscope (SEM). 

In this paper, we report on the properties of two films 
with different thicknesses, 150nm and lOOnm. The room 
temperature. X-ray diffraction pattern for the lOOnm 
film is displayed in Fig[TJa). Diffraction patterns for 
both films show similar features. A large signal is mea- 
sured for all relevant orders of the (£00) peak, confirming 
excellent monoclinic structure of VO2, found for r-cut 
sapphire substratesi^. No other vanadium oxide peaks 
were observed. The grain size distribution extracted 
from the SEM images (see (Figlljb)) and (c) was ~40nm 
for the lOOnm film, and ^130nm for the 150nm-thick 
film. Variation in grain size is due to the difference in 
sample thickness and probably also to intricate differ- 
ences in sample preparation, such as the exact deposition 
temperaturei^iiSii^i. 



B. Experimental Setup 

Thermal hysteresis properties of thin VO2 films were 
studied from analysis of resistance vs. temperature mea- 
surements by using a FORC procedure (see below) . The 
resistance characteristics were acquired in a home-made 
test bed by standard 4-probe measurements using a con- 
stant current source and a voltmeter at a sampling rate 
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FIG. 1: a) Room temperature. X-ray diffraction pattern of 
a lOOnm VO2 sample on a r-cut (1012) sapphire substrate. 
(^00) planes of the Ml phase are observed. Substrate peaks 
are marked along with instrumental background peaks (x- 
marked peaks), b) SEM micrograph (7 /xm side) of the lOOnm 
sample (small grain 40nm) and c) the 150nm sample (larger 
grain 130nm). 



of 10 Hz. A temperature controller was used to sweep 
the temperature at a constant ramp rate of 3K/min. In 
this experiment, it is crucial to avoid temperature oscilla- 
tions during sweeps, and especially during changes in the 
sweep direction. We used two Platinum temperature sen- 
sors to avoid any such effect: One sensor was positioned 
close to a resistive heater, and provided feedback to the 
control loop. The close vicinity to the heater reduces 
temperature oscillations in the sample which is farther 
away. The second sensor was in thermal contact with the 
sample. This sensor provided accurate measurements of 
temperature during the thermal cycles, along with the 
exact reversal temperature of each reversal curve, see be- 
low. 



C. FORC procedure 

A set of reversal curves for each sample was obtained as 
follows: starting at a temperature of 380 K (fully metal- 
lic state), the temperature was ramped down to a spe- 
cific reversal temperature, Tr. When Tr is reached, the 
sample is heated back to the fully metallic state. The 
path followed by the resistance, starting at Tr and end- 
ing at 380K {R{Tii,T)), is called a first-order reversal 
curve. Tn was varied between 370K and 260K in steps 
of 2K around the saturated regions; while in the mid- 
dle of the MIT, we used Tr steps of ~ 0.05K. Approxi- 
mately 200 FORCs were obtained for each measurement 
set. FORC measurements obtained with this procedure 
are called Heating FORCs, because temperature always 
increases for each Tr, see Figl^l^a). Similarly, we can 
define Cooling FORCs (Figl^^b)) by beginning with the 
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FIG. 2: A set of first-order reversal curves, a) If eating FORC. 
b) Cooling FORC. Solid lines are major loops, between fully 
insulating and fully metallic states. Doted lines show reversal 
curves log R [Tr, T) for three different Tr. 



fully insulating state, then heating up to various reversal 
temperatures, and finally obtaining a FORC for each Tr 
by cooling back to the insulating state. Herein, we show 
results for heating FORCs. To rule out any spurious ef- 
fects, due to sample degradation with measurements, we 
measured FORC sets with both increasing and decreasing 
Tr between sequential curves. No difference was found 
between the two procedures. For every set of data we 
calculated the FORC distribution, defined as the mixed 
second-order derivative^, 



p(Tr,T) 



2 dTRdT '■ 



(1) 



where r {Tr,T) = log (i? (Tr, T)). By calculating the 
FORC distribution for the log of the resistance there is 
an increase in sensitivity to features over the entire range 
of resistance, which is almost four dccades^S. The mixed 
second-order derivative eliminates the parts in the data 
set where resistance is constant under a change in T or 
in Tr. Thus, any non-zero value in the distribution cor- 
responds to irreversible parts in the hysteresis loop. 



III. RESULTS 

The heating FORC data sets measured for our VO2 
thin films are presented in Fig|3I^a) for the thin film sam- 
ple and FiglSt^d) for the thick film one. The major hys- 
teresis loops define boundaries of the reversal curves. The 
corresponding FORC distributions are depicted in the 
contour plot in Figl^^b) and[3Jc) for the thin and thick 
samples, accordingly. Figures[3jc) and[3l^f) arc 3D repre- 
sentations of these distributions, with identical z-scales, 
enabling direct comparison of magnitudes in the FORC 
distributions between the two samples. 

While the general features for both samples arc simi- 
lar, the FORC distribution contours enable us to easily 
compare various aspects of the two samples measured. 



and to also find some surprising features in the rever- 
sal curves, which are hard to identify otherwise. Three 
dotted lines are marked in Figl^Jb) and (e) to illustrate 
different reversal stages in the samples. The reversal tem- 
peratures of these curves are marked in Figures [3ja) and 
[31[d). These points represent the starting point of a spe- 
cific curve. 

The samples are insulating and fully reversible for tem- 
peratures below line 1 (in figures[3^b)-100nm sample and 
[3l^e)-150nm sample); hence, there is no reversal signal. 
This line corresponds to 327K in the lOOnm sample and 
324K in the 150nm one. Above line 1, a non-zero part 
evolves as Tr is increased, appearing as a tail composed 
of a positive-negative part, around T = 340K. 

Interestingly, at point 1 in the FORC curves (presented 
in figures [3Ka) and (d)), it seems there is no hysteresis 
for this reversal temperature. This feature is clear in 
the 150nm sample (Fig[3Jd)) where the hysteresis opens 
only 4 degrees higher, around 328K. But this point corre- 
sponds to the Tr, where the tail starts, i.e., where there 
is an irreversible feature in the FORC distribution. This 
feature and its meaning will be discussed in detail further 
below. 

As Tr increases, the tail joins a dip-peak structure 
which is the main part of the hysteresis, and includes 
most of the irreversible structure. The peak maximum is 
marked by line 2, around Tr = 334K for both samples. 
There are many quantitative properties of the hysteresis 
and reversal process that one can extract from the FORC 
procedure. It is clear that the peak for the lOOnm sample 
is wider and shifted to higher temperatures versus the 
150nm sample. Figure 3] illustrates this point, showing a 
cross-section that passes through the FORC peak and is 
parallel to T = Tr^ (see inset of FigU]). The square shaped 
plot is for the lOOnm sample and the star shaped one is 
for the 150nm sample. We find the maximum peak for the 
thin sample occurring at T=340.5K with magnitude of 
Pmax ~ 1.4x 10"^ and fuU width half maximum (FWHM) 
of 1.60. For the 150nm sample the transition is shifted to 
lower temperatures with the maximum at T = 334. 6K, 
a larger maximum of Pmax ~ 2.74 x lO"'^ and narrow 
distribution, FWHM of 1.30. The negative part of the 
FORC distribution in both samples has similar features 
regarding minimum and width. 

Above hue 3, in Figi^b) and E^d), the FORC dis- 
tribution is again zero, indicating that all R-T curves 
are completely reversible above this temperature and the 
sample is fully metallic. These temperatures correspond 
to 343K and 340K for the lOOnm and 150nm samples, 
respectively. 

Thus, the FORC analysis presented above provides 
much quantitative information of the samples studied. 
For example, Ref. reports that the hysteresis of a VO2 
thin film sharpens as the grains become larger, which 
agrees with our findings. By performing FORC measure- 
ments and analysis, one can quantitatively compare ad- 
ditional features of the hysteresis with various film char- 
acteristics. These include the center of the hysteresis at 
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FIG. 3: Complete set of First-order reversal curves, the corresponding contour plot distributions, and 3D view for the VO2 
films on r-cut sapphire substrate. Top row corresponds to the lOOnm-thick sample (40nm grain size) and the bottom row to 
f50nm-tick (130nm grain size). 
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FIG. 4: Profile lines parallel to T = Tr axis, passing through 
the pmax marked in the FORC distribution (inset) for the 
lOOnm sample (top) and the 150nm sample (bottom). Main 
panel - 100 nm sample (square-symbol) and 150 nm (star- 
symbol). Curves with same color map as the FORC plot. 



T and Tr, the FWHM. and temperature limits of the re- 
versible features. Such quantifying capabilities are valu- 
able tools for analysis and understanding of the relation- 
ship between these properties. 

More importantly, the FORC analysis enables us to 
observe features in the hysteresis, which are difficult to 
detect otherwise. As stated earlier, in both films there 
is a tail in the FORC distribution around T = 340K, 
extending to reversal temperatures lower than 325K. This 
feature is enhanced in the thicker (150nm) VO2 sample. 
Upon examining the major loops of the hysteresis curve 
in this sample, hysteresis begins to open between the 
major heating and cooling curves at ^ 329K. But, for 
reversal temperatures between 325K and 329K there is a 



non-zero value in the reversal curves around T = 340K, 
leading to a non-zero value in the FORC distribution; 
and only for Tu < 325A" the FORC distribution depicts 
a fully reversible region. 

This feature is similar to the one reported by Davis et 
al^ for magnetic hysteresis in Co/Pt multilayers. We 
follow a similar interpretation of our results. VO2 has 
been shown to have phase separation between metallic 
and insulating domains across the transitional^. If there 
are metallic domains which can persist to low tempera- 
tures, then when we measure a reversal curve (from these 
low temperatures) their effect on the resistance is negli- 
gible since they are small. But during the reversal curve, 
they may act as 'seeds' for the transition, and - through 
interaction with neighboring sites - they can modify the 
path of the MIT compared with the path taken from a 
fully metallic state. The change in path is manifested in 
the non-zero values and tail shape of the FORC distri- 
bution. 

We find further reinforcement for this scenario by mea- 
suring the resistance vs. temperature characteristics of a 
cooling curve on a 1 x 6 /xm^ VO2 junction with similar 
sample characteristics as the 150 nm sample, presented 
here. We have previously shown that this measurement 
enables us to identify the transition of a single domain 
from metallic to insulating^. Figure [S] depicts such events 
as clear jumps in the resistance even at 321. 6K. This 
means that down to these temperatures, metallic do- 
mains exist in the sample. The inset in FigIS] displays 
the full hysteresis loop of the device and the dotted line 
indicates the position of the jumps in the loop. It is clear 
that metallic domains persist down to temperatures out 
of the hysteresis loop. 

We wish to emphasize that the mere existence of metal- 
lic domains at lower temperatures should not change the 
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FIG. 5: Resistance vs. Temperature for a 1 x 6 firri^ VO2 
device. Main panel shows a zoom around 321K of the major 
loop presented in the inset (dotted line indicates 321K). Large 
jumps, indicating domain switching, appear at 321. 57K and 
321. 68K. 

shape of the FORC distribution unless these domains in- 
teract with the surrounding insulating domains, thus, af- 
fecting the course of the reversal curve. These interac- 
tions possibly play an important role through the entire 
phase transition and could be the driving force in the 
opening of the hysteresis in the VO2 MIT. In the follow- 
ing section, we derive such a phenomenological theory 
and show that it can capture the main features we see in 
the FORC measurements. 



IV. A MEAN-FIELD EFFECTIVE-MEDIUM 
THEORY OF HYSTERESIS IN VO2 

Before proceeding, let us briefly describe the more com- 
mon Preisach mode P°i'^^ . The Preisach model, in its sim- 



plest form, is based on the assumption that each grain 
has a transition temperature Tc±STc, where the plus and 
minus signs correspond to heating and cooling, respec- 
tively. Averaging over Tc and STc, the distributions of 
which serve as fitting parameters, one may obtain excel- 
lent fits between the model and the experimental data^i. 



In the model presented here, we look for the micro- 
scopic origin of the hysteresis loop in the inter-grain in- 
teraction, rather than assuming that it is already present 
for the single grain. The physical interpretation of our 
model is the following: each grain has its own critical 
temperature, given by its parameters (e.g. size, oxygen 
doping etc.). However, when the temperature exceeds 
the critical temperature of a specific grain, the grain still 
has to overcome an interaction-induced energy-barrier to 
change its phase. The energy barrier depends on the con- 
centration of same- type domains, which is the origin of 
the non-linearity required to generate hysteretic behav- 
ior. 



The calculation method is as follows: first we calcu- 
late the concentration p of metallic VO2 domains as a 
function of temperature in the presence of inter-grain in- 
teraction for both directions of the temperature scanning 
(this corresponds to the mean-field part of the calcula- 
tion). Once the concentration is calculated, it is used as 
input to calculate the resistance via the effective medium 
approximation (EMA)2S. The EMA relates the concen- 
tration of the different domains to the sample resistance, 
and has been successfully used to describe the MIT in 
VO2 filma^. The resistance can be evaluated by using 
the EMA via the formula^ 



R = 2 



1 + p (1 + 2) 2 + p (1 + z) - z / 4 (z - 2) RiRm + ii2+p{z-l)-z)Ri + {l+p{z- 1)) Rm) 



Ri 



Rm 



R^jR 



M 



where i?/, Rm are the resistances of an insulating and 
metallic grain, respectively, z is the number of nearest 
neighbors, and p, as stated above is the concentration of 
metallic domains. 

Let us start by calculating the concentration of metal- 
lic domains upon going from low to high temperatures- 
Pup (T')- Consider first a single grain (i) with its own 
critical temperature T^K If there were no interactions, 
then the probability of finding this grain in a metallic 

state would simply he pi = O (t — Tc^^ ^ , where O is 
the Heaviside step function. We now take into account 



the inter-grain interaction by assuming that the insulat- 
ing domains around the (i)-th grain generate a poten- 
tial barrier, which inhibits the transition. This means 
that in order for the grain to turn metallic, it is not 
enough that T > Tc , but there is also an energy 
barrier to overcome. If the average interaction energy 
is U, then the effective energy barrier under a mean 
field approximation is AE^ = (1 — pup)z[/, where p^p 
is calculated by averaging over all the domains (p^p)- 
The probability of escape from a potential barrier is 
proportional to the Boltzmann factor, and so we have 
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(l-Pup)zt/ 



In order to av- 



erage over all the domains, we assume that the critical 

(i) 

temperatures Tc of the different domains come from a 
Gaussian distribution averaged at the center of the tran- 
sition with a typical width AT. With this assumption, 
we find the equation for p^p , 



Pup = ^ ( 1 



exp 



(1 -Pnp) ZU 



(3) 

which is solved numerically. The cooling concentration 
Pdown may be calculated in the same way, now it is the 
fraction of metallic domains that generates the energy 
barrier. In order to obtain reversal curves, we calculate 
Pdown (T) from an insulating state at high temperature 
to Tfj, then use the concentration of metallic domains, 
and obtain pup (T^ Tr) the concentration of metallic do- 
mains for a heating FORC with a reversal temperature 
Tji. Once (T,T^) are obtained, the resistance is cal- 
culated via Eq. [51 

To fit the model to the experimental data, one needs 
the resistances of the metallic and insulating state. How- 
ever, these may be obtained experimentally from the low- 
and high-temperature data and are, thus, input parame- 
ters for the FORC analysis. The only fitting parameters 
are the temperature distribution width AT and the in- 
teraction energy U (assuming z = 4 for simplicity). Both 
of these parameters can be obtained only from the full 
hysteresis curve (i.e. solid line in Fig. 

In Fig. inja) we plot the theoretical curve R(T) for 
heating (blue line online) and cooling (red line online) 
measurements. The points are the experimental data. 
The parameters are i?/(T) = 36.6 exp (2||L5^ fl, Rm = 

25.7n, keATc = 2.26K and keU = 79.85K. Already, this 
curve shows good agreement between the data and the 
theory (noting that resistance changes in four orders of 
magnitude). In the inset of Fig. [^a) we plot the the- 
oretical FORC curve, which exhibits similarities to the 
experimental one, namely the "moon- like" shape and the 
peak-dip structure. 

In Fig. Ill^b-g) we plot the resistance curves for different 
values of Tf. The theoretical curve is the solid lines, and 
the points arc the experimental data. We emphasize that 
these curves, which exhibit good fit between the theory 
and experiment, were obtained without additional fitting 
parameters other then the two in Fig.[5Ja). 



V. SUMMARY 

In summary, we presented FORC measurements and 
analysis of the hysteresis in the temperature driven metal 
insulator transition of vanadium oxide thin films. We 
find that the FORC procedure, used extensively in char- 
acterizing magnetic systems is also a powerful tool for 
quantitative analysis of the MIT in VO2. 

The sensitivity of the FORC distribution to any ir- 
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FIG. 6; a) Resistance hysteresis curves obtained from tlie tlie- 
oretical (solid lines) and experimental data (points). The ar- 
row direction indicates heating and cooling. Inset: theoretical 
FORC curve, exhibiting similar features as the experimental 
one ("'moon-like'" shape, peak-dip structure), b-g) Resis- 
tance vs. temperature for different reversal temperatures Tr. 
Solid lines are the theoretical curves and points are experi- 
mental data. These curves were generated with no additional 
fitting parameters (see text). 



reversible feature enabled us to find that there are long- 
lasting metallic domains in the phase separated VO2 film, 
existing down to temperatures where the sample appears 
fully insulating. Moreover, we find a signature in the 
FORC analysis, indicating that these domains interact 
with the surrounding medium, modifying the path of the 
transition. 

In this spirit, we developed a mean field theory to cap- 
ture the hysteresis. Here, the driving force in charge 
of opening the hysteresis of the MIT is the interaction 
between nearest neighbor domains and the two fitting 
parameters having physical meanings. This is enough to 
reconstruct the main experimental features. 

We hope this paper will encourage the use of the FORC 
method in other systems where appropriate, and that it 
will add to the theoretical discussion regarding the origin 
of hysteresis in general. 
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